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ABSTRACT 

With the advancements in engineering and technology and with the advent ofthe fourth industrial revolution 
made the solution for complex and dynamic problems easier. Many researchers had incredibly investigated about the 
problems, which are outcomes ofprosthetic leg, attached to a human being especially at working conditions and still the 
problem pertains till date. This might be one of the crisp reasons which insisted& necessitated to undergo this project. 
In this project, considering a set of persons aging 19-28 were made to walk on a tread mill. It is a tedious process to take 
the hip and knee angles physically using various sensors and its very complicated process to generate torque values and 
finding the optimal torque values for different persons converting physical variables to digital variables. To overcome 
this problem advanced software’s like Kinovea is employed for measurement of angular cycles for gait analysis and 
based on these values machine learning algorithm is employed to find the accuracy and better prediction of torques for 
different disabled persons. 

KEYWORDS: Frames per second (fps), Gait Analysis, Multivariable Linear Regression (MLR) & Root Mean Square 
Error (RMS) 
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1. INTRODUCTION 

In early 1990’s it was very difficult to measure the hip angle as well as the knee angle of a specially abled 
person because there it needed a physical mounting of sensor on the amputee to calculate the knee angle as well as 
hip angle for different speeds and different walking cycles [1][2]. Though the employability of dart fish sensors as 
well as required values, still there is an uncertainty about mounting of sensors and calibration of sensors in 
measuring the angular ankle cycles. Gait analysis for these values for an amputee is not that accurate and 
ergonomical, these constraints had made a huge necessity to go far advanced software’s like kinovea which works 
on video motion analysis [3] [4]. Kinovea is best among the present available motion analysis software’s for 
capturing real time values at every instant of time and it is available in open source [5]. Generally in the case of 
developing nations like India and many other countries having the huge demand of exoskeleton for specially 
disabled amputees for various reasons like accidents and genetic problems, based on current population and 
increasing number of exoskeletons in the field of transfemoral prosthetics. It necessitated for cost optimization as 
well as better efficient and accurate test result values from various sources. Employing automation that is highly 
accurate and strong, efficient in programming and deriving outputs makes a difference. With the rapid growth of 
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industrial revolutions the latest software’s like python is used for reducing the complexity of problem, by easing the 
mathematical empirical calculations for measuring maximum cycles. Later on including the methodologies of machine 
learning procedures gives better accuracy and efficient results for designing and developing exoskeletons in more 
ergonomical & efficient manner. 

In present world populating more than 7 billion persons around the globe, more than 17 to 20% in a billion of 
them were in the need of prosthetics and exoskeletons for various reasons like aging, amputing because of accidents, 
severe illness and physical stress. For procuring gait analysis needs calculation of torque values by using Lagrangian Euler 
formulation, one of the best methods of analysing and computing the complex structure and the behaviour of serial link 
chain which is efficiently done by Lagrangian equation, this method provides easier way of solution. 

2. METHODOLOGY 

This paper depicts mainly on the prediction of output torques for the gait cycle that has to be performed to the 
amputee using machine learning techniques. In this paper, the torque governing equation is derived from the Lagrangian 
mathematical model. Various independent inputs like height, mass, weight of the amputee is obtained through video 
motion analysis. Acquired data is processed using advanced computational software and the prediction of torques can be 
obtained by employing machine learning algorithms. This has been pictorially represented with the flow chart that 
describes the series of steps that has to be performed to obtain the output torques for amputees. 
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Figure 1: Flow Chart for Prediction of Torques 


2.1 Lagrangian Model 


Lagrangian equation is mainly based on generalized forces, as well generalized coordinates are defined as the 
angles and distances considered for revolute, prismatic joints [6]. The torques or forces that were associated with the joint 
actuators were calculated using generalized forces. In general terms for a physical system the total number of generalized 
coordinates are considered, calculated depending on the position of the mass of degrees of freedom. This modelling 
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technique is based on the lumped parameter model with internal components of point masses which are located at centre of 
gravity of the generalized links. The Lagrangian function is represented with ‘L’ and is defined as the difference between 
kinetic energy (K) and potential energy (P) of a system. 

L= K-P (2.1) 

The Lagrangian Euler formulation is described with a set of system variables. For the generalized displacement 
‘q’ for a joint variable the Lagrangian Euler formulation is obtained by 

— (-=x,; fori= 1, 2,-n (2.2) 

dt v dqr dqi 15 ’ ’ v 7 

For the condition Tj =0, states that joint T is static whereasTi ^O, the system is dynamic and is modified at various 
instance of time. 

For a human leg, the degrees of freedom are mainly dependent on two functional parameters i.e knee joint and hip 
joint. Hence the system is of 2 DOF and it requires two torque governing equations to be calculated. 



Figure 2: Lagrangian Model for 2-DOF System 


The joint variables were represented as 0i, 0 2 and link lengths were represented as Li and L 2 mass of the links 
were represented as ml and m2 respectively. The linear velocities were represented as Vi and v 2> angular velocities are 
represented0i,0 2 respectively. For link 1, substituting and generalizing all the process variables the kinetic energy of link 1 
is obtained as 

Ki=imivi 2 +4iC0i 2 =i miLi 2 0i 2 + ^ miLi 2 0i 2 = \ m x U% 2 (2.3) 

The potential energy for link (1) is given as 

Pi = \ migL^sine! (2.4) 

Considering link (2), the Cartesian coordinate position values(x 2 , y 2 ) are given as 

x 2 = Licos0i + ^L 2 cos (0i+0 2 ) v (2.5) 
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y 2 = LisinQi + ^L 2 sin (0i+0 2 ) (2.6) 

Differentiating those equations and simplifying the magnitudes of velocities of link 2, the final kinematic equation 
is obtained as 


K 2 - -m 2 v 2 + -1 2 0)2 2 (2.7) 

The potential energy for link 2 is obtained as 

P 2 =m 2 gL 1 Si+ - m 2 gL 2 Si 2 (2.8) 

By substituting the kinematic as well as potential energy the torque equation for link (1) is derived as: 


^i-[(“ mi+m 2 )Li 2 +- m 2 L 2 2 +m2LiL 2 C 2 ]0i+m 2 [-L 2 2 +-LiL 2 C 2 ]0 2 -m 2 LiL 2 S 2 0i0 2 - - m 2 LiL 2 S 2 0 2 2 +(- mi+m 2 )gLlCl+- 
m 2 gL 2 Ci 2 (2.9) 

Similarly for link (2) the governing torque equation is derived as 

l 2 = m 2 [i L 2 2 + \ L1L2C2]0! + i m2L 2 2 0 2 + i m 2 L 1 L 2 S 2 9 1 2 + \ m 2 gL 2 C 12 (2.10) 

From (2.9), (2.10) governing equations the empirical relation for finding the torque values for a 2 DOF prosthetic 
system is given as: 



( mJKl+Ci) m JKl + - z C 2 )\ 

01 

+ 

—m x l\S 2 Q \0 2 — ^n\l\S 2 ^ 2 

\m 2 l 2 2 {\ + \c 2 ) \m 2 l 2 ) 

. 02 . 


-m 2 l 2 S 2 Qi 


\ m idliCi2 + 

2 m 2dhCl2 


2.2 Kinovea 


( 2 . 11 ) 


Mounting of physical sensors on the amputees for obtaining gait analysis and finding torque values is never been 
that easy. To overcome this problem, with the advancements in the field of engineering and technology simulation 
software’s were employed, which performs objectification of human movements. Kinovea is the best open source software 
which is a low cost technology that enables motion analysis for 2D kinematics parameters. Kinovea software enables the 
user to analyze the movement of any subject, calculation of displacement and mainly employed in calculating gait analysis 
[7]. Kinovea enables the analysis of angles and spatiotemporal parameters. It enables the user to measure angles and 
distances with respect to frames. It has high reliability in measuring hip, knee and ankle measurement values with different 
velocities with respect to time. Kinovea employs point tracking function which enables to obtain data without any training 
of the data to system, which is reliable and efficient. Kinovea is easy to use as well as portable and highly accurate in 
obtaining kinematics data in real time scenario. It employs image capturing procedures, frames calibration, image 
digitalization and finally exporting of data to spread sheet by data extraction and transformation techniques for statistical 
analysis of gait cycle. In this research kinovea 0.8.15 version, which is an open source downloadable format for any 
operating system is used. 
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Machine learning is the process where a computer can learn themselves from the data provided as input by the 
process of continuous training and testing and making the system intelligent[8]. With the help of machine learning 
concepts, there is no necessity of explicit programming for the huge amount of data and the tasks that are to be 
programmed in the earlier cases [9]. The performance of machine learning is measured in terms of accuracy of the 
algorithms in the process. 
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Figure 3: Machine Learning Block Diagram 

The learning concepts in machine learning are broadly classified as 4 types. They are: 

• Supervised learning. 

• Un-supervised learning. 

• Semi-supervised learning. 

• Reinforcement learning. 

In this paper by employing supervised learning process the required output data is generated. According to the 
inputs the best suitable algorithm for processing the data for testing and training conditions is multiple linear regression. 
For programming in machine learning, opted python script, as python is one of the best open source software which is 
available in market and is having high performance and strong programmability of machine learning concepts. Employing 
spyder, cross platform integrated development environment for scientific programming language python which is open 
source IDE tool, available in anaconda navigator app. By downloading anaconda navigator app and launching the spyder 
software, programmed the given sets of input training data using multiple linear regression algorithm. 

3. DISCUSSION & COMPUTATIONAL ANALYSIS 

Experimental Setup 

The data has been collected from 50 healthy student’s age ranging from 19-28 years of age [10]. The persons were 
not prone or subjected to any injuries and abnormalities. The data was collected from the persons during a study state 
controlled speed at 2 km/hr on a treadmill using video motion analysis using kinovea[l 1]. The lens opted for coupling the 
data during the steady state controlled walking speed is Sony F57 E-mount 28-135mm F4 cineservo 180frames/sec(fps) 
with the shutter speed 1/200 s, distanced at 3.5mt from the treadmill and recorded the walking stances of the person. Three 
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colour node indicators were attached to the hip, knee and ankle to calculate the respective angle measurements. The 
recorded data is for the duration of one minute and fps is at the rate of 60 fps. 



Colour node 
indicators 


Treadrnill 


Figure 4: Experimental Setup on a Treadmill 


3.2 Data Retrieving 

The captured data is calibrated for one complete cycle which includes gait cycle i.e. stance phase and swing phase 
and recorded the values with reference to the colour strips perpendicular to the sagittal plane axis and tabulated the values 
hip, knee with respect to the angle [12]. This process is carried out by loading the video content to kinovea software by 
which motion analysis happens and accurate results were obtained. Basically the uploaded video content for duration 60 
seconds is calibrated to the time taken for a person to complete one gait cycle. This procedure is repeated for 3 cycles from 
the same subject for the gait analysis by considering 78 fps for one gait cycle and tabulated hip, knee angles [13]. The same 
sequence of steps has been performed for 50 different subjects of different age groups heights weights and noted down the 
respective values. 



(C) (d) 

Figure 5: Hip & Knee Angle Measurement in Kinovea (A) Heel Strike, 
(B) Toe Off, (C) Mid Swing, (D) Heel Strike 
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(c) (d) 

Figure 6: Linear Velocity Measurement in Kinovea (a) Heel Strike, 
(b) Midstance, (c) Toe Off, (d) Mid Swing 


That recorded values has been grafted from spread sheet with provision provided by the kinovea software, which 
made tabulation of values simpler. The data which is retrieved from kinovea software is the source for obtaining torque 
values from the Lagrangian Euler model. The data obtained through kinovea are 

3.2.1 Hip & Knee Angle (Rad) 

Hip and knee angle data is obtained through kinovea by direct subject measurement, i.e., the reference point is 
taken from the strip attached to the hip of the subject (person) with respect to sagittal plane at different instances of time, 
by this hip angle(0i) data is obtained similarly the knee angle data is obtained by taking the reference point from the strip 
attached on the knee of the subject (person) with respect to sagittal plane at different instances of time, by this knee 
angle(0 2 ) data is obtained and recorded. 

3.2.2 Hip & Knee Angular Velocity (Rad/Sec) 

Hip angular velocity (0i) data is obtained by differentiating the values of hip angle (0i) with respect to time. Knee 
angular velocity (0 2 ) is processed by differentiating the data of knee angle (0 2 ) with respect to time. 

3.2.3 Hip & Knee Angular Acceleration (Rad/Sec 2 ) 

Hip angular acceleration (6?i) is calculated by differentiating the values of hip angular velocity (0i) with respect to 
time. Similarly knee angular acceleration ( 0 2 ) is obtained by differentiating the data of knee angular velocity (0 2 ) with 
respect to time. 

The data tabulated from above values has been normalised and represented in the form of graphical canonical 
form below 
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Figure 7: Depiction of Hip Angle, Hip Angular Velocity, Hip Angular Acceleration 



Figure 8: Depiction of Knee Angle, Knee Angular Velocity, Knee Angular Acceleration 
3.3 Data Processing 

The data obtained from kinovea software through spreadsheets is been tabulated for measurement of hip and knee 
angles of different aged persons of different mass and heights. The calculation procedure for obtaining the torque of the 
motor that has to be employed in the prosthetic is very complex and time consuming as well as error prone [14]. The 
calculation procedure of final torque involves the substitution of values obtained from kinovea software are hip angle, knee 
angle, linear velocity at hip and knee with respect to time. Along with these values, mass of thigh (mO, mass of shank (m 2 ), 
length of link li(thigh), length of link 1 2 (shank) and 0i ,0 2 , 

0i,0 2 ,0*i, 9 2 are calculated 

Where 0i ,0 2 = angular displacement of first and second link 

0i,0 2 = angular velocity of first and second link 

0i, 6 2 = angular acceleration of first and second link 


0 12 = 0j + 0 2 , Si=sin0i, S 2 =sin0 2 , Si 2 =sin0i 2 , Ci=cos0i, C 2 =cos0 2 and Ci 2 =cos0i 2 
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This whole procedure is cumbersome, complex and there is a chance of occurring of error while calculating the 
final torques. 


This procedure has to be perform for 50 persons for different parameters at the rate of 72 values. This procedure 
has to be repeated for 3 cycles. So in order to overcome this, employed the strongest simulation tool which can perform 
this whole complex structure and gives accurate results, that software is python, available in open source cross platform 
IDE and calculated the output torque values. Employing special packages in python like numpy and math libraries eased 
finding the output torque values. 


3.4 Multiple Regression 

The main significance of multiple regression algorithm is to get the relationship between inputs and output where 
the inputs are having more number of independent variables and only one dependent output variable. For this research, 
having independent variables like age, heights and weights which are multiple in numbers and the output dependent 
variable, what has to be obtained is torque. As it was nearly impossible to calculate the output torque with respective 
multiple number of independent input variables, employed multiple linear regression technique and formulated the output 
torque with respect to input values [15] [16]. 

Considering the multiple linear regression model having independent predictor variable ‘X’ X=(x ls x 2 , x 3 ... x n ) 
and dependent output variable ‘y’ is given as 


y = Po+PiXi+...p x x k +e 


(3.1) 


For ‘i’ number of iterations y- x is given by 


yi = p 0 +PiXii+...+(3 x x ik + Ciwhere i= l,2,3...n 


(3.2) 


The sum of residual square errors is given by 


X e ' 2= X (yi-Po- X Pi X ij ) 2 


(3.3) 


That is more simplified as 


X £i 2 = e'e ^(y-XpPly-XP) (3.4) 

i =1 


The multiple linear regression technique employed the importing of numpy and pandas packages which are inbuilt 
and available in python source module spyder from anaconda app, by downloading anaconda navigator and installing 
spyder software imported the data in the acceptable format for data processing through numpy and pandas. The main 
importance of numpy is to calculate and predict the output values from the huge data sets that data has been imported in the 
form of data frames for processing in pandas. This has a very significant and promising role in machine learning and deep 
learning techniques. Employing a special graphic library i.e. matplotlib helps in plotting the values after processing of the 
dataset shown in figures (9),(10). A special package namely scikit-learn is employed for performing pre-processing tasks 
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like transformation of input data, standardizing of data for regression and classification operations [17] [18]. By using 
python script prediction operation using multiple linear regression algorithm is done and torque values were calculated as 
well as predicted from the input data. 
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Figure 9: Scatter Plot for Age, Height, Weight and Hip-Torque Weight 
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Figure 10: Scatter Plot for Age, Height, and Knee Torque 

Accuracy 

The predicted output value from multiple linear regression has to be validated by the standard means of 
calculating the error [19]. For this root mean square error method to be performed has a unique reason, imbibing RMS 
method, the accuracy of the machine learning algorithm to measure the difference between the predicted values and the 
actual values can be obtained using the equation (3.5), calculation of rms value for the output test and predicted cycles is 
tabulated. That rms value has been compared with the output predicted value and the error rate is ±4.069 for the predicted 
hip torque and ±1.496 for the knee torque shown in figures (11), (12). 
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RMS = Jig (y r y) 2 (3.5) 

n = number of test data elements (01, 2 ...n) 

Yi = Actual value 
y = Predicted value 
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DataFrame 
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y_pred 
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(29,) 

[351.5146865 358.66826417 363,38816612 ... 352.14526877 
345.36661367 ... 

y_test 

float64 

(29,) 

[346.1433163 359.7943807 373.8524167 ... 356.3762925 348.6... 
368. ... 

y train 
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(43,) 

[379.4663331 313.3116432 322.6175127 ... 368,7366468 425.3... 


Figure 11: Rms Error Value for Hip Torque Prediction 
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,, 26,25491534 27.37... 
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— 

,, 26.46668286 29.47... 


Figure 12: Rms Error Value for Knee Torque Prediction 


4. RESULTS & CONCLUSIONS 

The results obtained through machine learning technique by employing the multiple linear regression algorithm is 
been recorded. The predicted output torque values obtained after processing of data sets by python script through spyder 
software is shown in figure below. 
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Figure 13: Predicted Hip Torque 
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-XZ-fi- 
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Figure 14: Predicted Knee Torque 


New data set with three different values of age, height, weight of a person respectively were fetched, the data that 
is considered for the prediction isn't been recorded in the tests and trained data set model. With the help of multiple linear 
regression algorithm the predicted output torque for the given is of having tolerance ±4.069 and +1.496 values of the 
original calculated value. The data obtained from kinovea is compared with the previous methods of data retrieving 
processes. The data when compared is more or less similar to the previous methods and hence the processing of that data 
by machine learning algorithm and the prediction of the torques were said to be accurate and hence this method is 
recommended for prediction and calculation of torque values through gait analysis. This paper shows that the output 
predicted torque values were almost equivalent to theoretical calculations. This method reduces the complexity of finding 
torque values from the gait cycle and reduced the redundancy of output values this implies that by formulating and taking 
the same procedural steps can be used for obtaining torque values for a person who needs transfemoral gait analysis and 
the amputees for better output processing of torque which indeed helps in selection of appropriate motor in prosthetics. 

5. FUTURE SCOPE 

With the advancements in science and technology by industrial revolution 4.0, there could be development in 
much more efficient video motion capturing analysis software's, which intend helps in data capturing and retrieving 
process. The emerging growth and development in the fields of deep learning concepts could increase the development and 
enhancement of better algorithms, optimizing functions for prediction of values with minimum average loss in testing and 
training data sets. This helps a lot in mere future for better prediction and generating best outputs for the given data. All 
these techniques put together helps in determination of torque values in better way by employing data science concepts. 
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